Constraining the extra heating of the Diffuse Ionized Gas in the 
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ABSTRACT 

The detailed observations of the diffuse ionized gas through the emission lines Ha, [Nil], and 
[SII] in the Perseus Arm of our Galaxy by the Wisconsin Ha Mapper (WHAM)-survey challenge 
photoionization models. They have to explain the observed rise in the line ratios [Nil] /Ha 
and [SII] /Ha. The models described here consider for the first time the detailed observational 
geometry toward the Perseus Arm. The models address the vertical variation of the line ratios 
up to height of 2 kpc above the midplane. The rising trends of the line ratios are matched. The 
increase in the line ratios is reflected in a rise of the temperature of the gas layer. This is due 
to the progressive hardening of the radiation going through the gas. However an extra heating 
above photoionization is needed to explain the absolute values. Two different extra heating rates 
are investigated which are proportional to n° and n 1 . The models show that a combination of 
both are best to explain the data, where the extra heating independent of density is dominant 
for z > 0.8 kpc. 

Subject headings: radiative transfer, methods: numerical, Galaxy: halo, ISM: structure 



1. Introduction 

Containing typically half the mass of ionized 
hydrogen in galaxies, the Diffuse Ionized Gas 
(DIG) is visible as an extended Ha emitting layer 
in our Galaxy (e.g. the WHAM-survey) and 
in many other galaxies (see e.g. Tullmann & 
Dettmar (2000), Collins & Rand (2001), Otte et 
al. (2001), Hoopes & Walterbos (2003)). Stud- 
ies of emission line ratios such as [Nil] /Ha and 
[SII] /Ha provide information about the physical 
conditions of the gas. Simple energy estimations 
(Reynolds (1990), Reynolds (1993)) favor O and 
early B stars to be responsible for most of the DIG. 
3D models using various methods (e.g. Miller & 
Cox (1993), Dove & Shull (1994), Wood & Loeb 
(2000), Ciardi et al. (2002), Wood et al. (2004)) 
showed that it is possible for ionizing photons from 
O stars to penetrate from the midplane into the 
halo. Wood & Mathis (2004) noted that the line 
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ratios increase with distance from the midplane 
due to the progressive hardening of the radiation. 
So far photoionization models made no specific 
attempt to model the trends of the line ratios. 
Models by Mathis (2000), Domgorgen & Mathis 
(1994), Sembach et al. (2000), Bland-Hawthorn 
et al. (1997) used volume average models to ex- 
plain the observed data. The analytical approach 
by Haffner et al. (1999), referred to as Haffner99 
in the rest of this paper, treated the dependence of 
the line ratios with height. Haffner99 and its fur- 
ther application to other galaxies by e.g. Collins 
& Rand (2001), Otte et al. (2001), and Miller & 
Veilleux (2003) gave evidence that an additional 
heating source is needed in order to explain the 
rise of the line ratios with increasing distance z 
from the midplane. 

We are constructing photoionization models in 
order to examine the trends in the observed line 
ratios and if photoionization can heat up the gas 
sufficiently in order to explain the data. We in- 
troduce specific extra heating terms, extra means 
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in addition to photoionization, and discuss their 
properties. In the following section 2 we intro- 
duce the observations of the Perseus Arm to which 
our models are compared. Section 3 deals with 
the model parameters and discusses geometry and 
sight line effects taken into account. The models 
are compared with the data in the next section. 
The last section summarizes the results. 

2. Data of the Perseus Arm 

We are using the data taken from the Wiscon- 
sin Ha Mapper (WHAM) (e.g. Haffner99, Haffner 
et al. (2003)) survey which were kindly pro- 
vided by Ron Reynolds and Matt Haffner. The 
WHAM-survey mapped the northern sky in Ha 
with declinations of 5 > —30°. The Perseus Arm 
(-35° < 5 < -11° and 120° < I < 150°) was ad- 
ditionally mapped in [NII]A6583 and [SII]A6716. 
At each pointing an averaged spectrum with a 
beam of 1° is measured with a velocity resolu- 
tion of 12 km s _1 . The emission of the Perseus 
Arm can be separated in velocity space from the 
local emission (Haffncr99) for galactic longitudes 
120° < 1 < 150°. This was performed by in- 
tegrating the line emission in the velocity range 
— 100 km s _1 < v < —20 km s _1 , no line fitting 
was performed. We are using the intensity of the 
Ha line as well as the line ratios [Nil] /Ha and 
[SII]/Ha. The sensitivity limit of 0.1 Rayleighs 1 
results in an observed vertical height of up to \z\ 
= 2 kpc assuming a distance to the arm of 2.5 kpc. 

3. Model parameters 

We use the spectral simulation code CLOUDY, 
version 96.00 (described by Ferland (2002), (2000), 
(1998)) to model the DIG. CLOUDY determines 
the physical conditions by balancing the heating 
and cooling rates, so that the energy is conserved. 
The results of the models are compared to the ob- 
served emission line ratios [Nil] /Ha and [SII]/Ha, 
and the gas temperature as derived from [Nil] /Ha. 

In order to realize a model describing the DIG 
certain parameters have to be specified: The ion- 
izing spectrum of the source is a composition of 
three different stellar temperature: 56% from T = 
35000 K, 12% from T = 40000 K, and 32% from 
T = 45000 K, as used in Mathis (2000) and Wood 

'lR= 10 6 /4tt photons cm~ 2 s" 1 sr" 1 



ct al. (2004). The WMbasic models (Pauldrach 
ct al. (2001)), which include N-LTE effects, X- 
ray emission from shocks within stellar winds, are 
used as the ionizing spectra. The luminosity of the 
source is chosen in such a way that the observed 
run of the Ha intensity is matched, as shown in 
Figure 1. 

The density structure is exponential, as derived 
from the observed Ha intensity, with a scale height 
of 1 kpc and a midplane density of 0.2 cm~ 3 . The 
density is in clumps with a filling factor of 20%, i.e. 
only 1/5 of the volume is filled with this plasma. 
The geometry of the ionized gas is chosen to be a 
plane parallel layer, i.e. the ratio of the depth of 
the cloud to the distance of the illuminated face 
to the ionizing source is smaller than 1/10. 

We are matching the observed Ha intensity of 
the observations (Figure 1) by placing the illumi- 
nated face of the cloud at a z-height of 1 scale 
height of the Lockman-layer (300pc) assuming a 
density law of n = 0.1 cxp(— z/0.3) as in Miller & 
Cox (1993). This is done after the model is calcu- 
lated as otherwise the condition of plane parallel 
illumination of the cloud cannot be fulfilled. The 
information of the actual position of the ionizing 
stars are effectively removed, consistent with the 
picture of having the DIG being ionized by radi- 
ation leaking out of the Lockman-layer. The in- 
tensity gradient of Ha is matched as well as the 
estimate of the hydrogen ionizing photon flux by 
Reynolds (1990)): 4>dig > 5 x 10 6 hydrogen ion- 
izing photons cm~ 2 s . In order to compare the 
models on a common basis, which is important 
as extra heating is dominating for large z-heights, 
we choose to let the ionization structure be the 
same for all models (see Figure 5). This is in 
accordance with the idea that the extra heating 
affects only the temperature of the gas. The for- 
bidden levels of nitrogen and sulphur can then be 
more easily excited by collisions with the electrons 
which in turn elevates the line ratios [Nil] /Ha and 
[SII]/Ha. Figure 1 shows that the intensity gradi- 
ent of Ha varies only a little for the different mod- 
els. The scale heights for the models are slightly 
different (<; 10%) for the models with an extra 
heating, which can be explained by assuming that 
the DIG-layer is in pressure equilibrium. This ef- 
fect was also noted by Wood & Mathis (2004). 
The ionization parameters (U) of all models lie in 
the narrow range between log(U) = -3.0 and -3.1. 
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The ISM composition of CLOUDY is used with 
N/H and S/H set to the values in Haffner99 (N/H 
= 7.5 x 1CT 5 , S/H = 1.86 x 1(T 5 ). Graphite and 
silicate grains with the size distribution used for 
the ISM (see HAZY, Ferland (2002)) are present 
in the gas and account for less than 10 % of the 
global heating. The inclusion of PAHs give only 
variations < 3% for the line ratios and have only 
a small effect on the heating balance. The inter- 
action with cosmic rays is taken into account as 
described in Ferland & Mushotsky (1984). 

A crucial factor is the consideration of the line 
of sight. The observations give information about 
the line ratios at different positions orthogonal to 
the source of the ionizing radiation. The correct 
local line ratios are represented by the ratios of the 
volume emission coefficients ey: 4 1 = I ^ eVl « 

h J dz e V2 

— . An important issue is the treatment of the 

e V2 * 

line of sight to the Perseus Arm due to our posi- 
tion in the Milky Way as shown in Fig. 2. Each 
pointing of the WHAM survey contains contribu- 
tions from different \z\ heights above the midplane 
(e(zl) to e(z2)), this effect is taken into account 
by integrating over the particular sight line in the 
models. The distance to the Perseus Arm is as- 
sumed to be 2.5 kpc and the thickness of the arm 
to be 1 kpc as quoted in Haffner99. Moreover, the 
observations are the average of a beam size of 1° 
which means that we have to take into account 
an additional integration over different z heights, 
on average 45 pc. These effects of geometry and 
observational smearing are taken into considera- 
tion here for the first time. Fig 3 shows the ef- 
fect of the beam smearing and line of sight ge- 
ometry on the calculated line ratios for z-heights 
up to 2.5 kpc, which have to be considered for 
the observed heights of 1.8 kpc. These effect in- 
crease the 'corrected' line ratios by at most 15% 
and their trend is altered for high z -heights. The 
'corrected' model without an extra heating seems 
to suggest that [Nil] /Ha is not increasing after 
about 1.6 kpc, but the 'uncorrected' model shows 
a further increase due to the hardening of the ra- 
diation. The result of the gas temperature (sec- 
tion 3.2, Fig. 4) is even more influenced by this 
issue. The discussion of the line ratios therefore 
need to include an appropriate treatment of the 
line of sight geometry and the beam smearing to 



take these effects into account. 

3.1. Extra heating 

Additional heating to photoionization is in- 
cluded which is proportional to n 1 and n° in 
accordance with Reynolds et al. (1999), their 
factors Gi and Gi respectively. We choose rates 
in the same range as in their paper: G\ = 1 x 
10 _25 ergcm~ 3 s _1 and G2 = 5x 10 _27 erg cm~ 3 s _1 . 
The heating-cooling balance can then be written 
as either G + G\/n e = A or G + G 2 /n 2 = A. 
The heating due to photoionization is given by 
GquI and the cooling by An 2 . The inclusion of 
an extra heating source rises the gas temperature 
of the models, at the same time the ionization 
structure varies only slightly. The temperature of 
the CLOUDY models are calculated by heating- 
cooling balance. The extra heating therefore in- 
creases the temperature, having more pronounced 
effects at larger z heights as the photoionization 
heating rate decreases like n 2 as shown in Fig. 4. 
The graph is explained in more detail in the next 
section. 

In our models the ionization structure is nearly 
unaffected by the inclusion of an extra heating 
source, sulphur is slightly more effected than ni- 
trogen. This behavior is the basic assumption 
for the extra heating in Reynolds et al. (1999). 
However they are assuming a constant ratio of 
N+/N, whereas the models show a dependence 
on z height, this is expected as the radiation gets 
progressively absorbed. Figure 5 shows ionization 
structure and the change dependent on the differ- 
ent heating rates. Hydrogen is nearly fully ionized 
throughout all models which is a basic character- 
istic of the DIG. 

3.2. Comparison with Observations 

In Fig. 6 and Fig. 7 the models differ by the 
type of extra heating, ranging between models 
without extra heating and a rate of G\ = 1 x 
10 _25 ergcm -3 s _1 and G2 = 5x 10 _27 erg cm~ 3 s _1 . 
There was no fitting done in order to match the 
models with the observations. No 'best-fit' ex- 
ists and therefore the models are independent of 
the observational quality, individual spectral fea- 
tures, or small scale variations which cannot be 
reproduced with a smooth density distribution. 
As the data for small \z\ heights are contaminated 
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with radiation from the midplane and dust ab- 
sorption it is convenient to consider line ratios for 
|z|-heights above 0.8 kpc to be 'pure' DIG. This 
is also the range for which the Ha scale height 
was determined. The models show lower values of 
[Nil] /Ha and [SII]/Ha for the lower z heights as 
doubly ionized nitrogen contributes 35% and dou- 
bly ionized sulphur even 80%, as seen in Fig. 5. As 
a consequence [Nil] and [SII] are weaker. Figure 6 
shows the development of the line ratios with |z| 
height, both the data and the models show an 
increase with \z\. The trend in the line ratios 
is matched by the models even without an extra 
heating source. 

The [Nil] /Ha line ratio above 1 kpc can be ex- 
plained with the models including an extra heating 
rate. The modeled [SII] /Ha ratio is up to a factor 
of two below the observations, however the theo- 
retical uncertainty for the [SII] line is very high 
as dielectronic recombination is an important pro- 
cess in the DIG. As the corresponding recombina- 
tion coefficients are not known (see discussion in 
Ferland et al. (1998)), the results of the models 
have to be handled with care. Our models use 
the KLUDGE approximation (Ferland (2002)). 
Models without dielectronic recombination have 
[SII] /Ha decreased by 50%. If the rate is doubled 
then [SII] /Ha is increased by 50%. 

Figure 4 shows the gas temperature of the mod- 
els, the increase of the line ratios with z height is 
due to the progressive hardening of the radiation 
as the photons go through the gas layer. In or- 
der to explain the observed line ratios an extra 
heating source is however needed which does not 
alter the general shape of the predicted line ra- 
tios but elevates the line ratios. As [Nil] and 
[SII] arc forbidden lines which get collisional ex- 
cited by electrons, an increase in gas temperature 
increases the amount of electrons capable to ex- 
cite the singly ionized nitrogen and sulphur ions 
which can then decay by emitting the emission 
lines in questions. The gas temperature is de- 
ducted from the observations through the rela- 
tion - following Reynolds ct al. (1999): = 

1.84 x 10 5 (^) (^)~V 39 exp(-2.18/T 4 ). 
Assuming that N + /N = H + / H and the abun- 
dance as stated in section 3 gives J ^" ] = 
13.75 T 4 - 39 exp(-2.18/r 4 ), T 4 = T/IO 4 ^." We 
are not using the collision strength of singly 



ionized nitrogen from Reynolds et al. (1999): 
N+ : n D) = 2.28 • T 4 026 , Aller (1984), 
but the data from Stafford et al. (1994): N+ : 
ft fP, 1 D) = 3.02 • T 4 ~ 001 . This leads to the dif- 
ferent coefficients and temperatures on average 
250 K below the values of Reynolds et al. (1999). 
For singly ionized sulfur we use the data from 
Lanzafame et al. (1993) instead of Aller (1984) 
as in Reynolds et al. (1999). The plot also shows 
the impact of the line of sight geometry and beam 
smearing as well as the consideration of the ioniza- 
tion structure of the models. As the line of sight 
leads to higher [Nil] /Ha ratios (up to 15%,) this 
effects, when accounted for, lower temperatures. 
The ionization structure, i.e. N + /N w H + /H is 
only valid for z > 800 pc, leads to higher estimates 
of the temperature. The combined effect is the el- 
evation of the derived temperature for z < 800 
pc and above that lower values by about 200 K. 
The models without an extra heating source have 
temperatures too low, an extra heating source in- 
dependent of density shows very good agreement 
with the observed temperatures. The extra heat- 
ing oc n 1 seems to best fit the data up to 800 pc, 
then the extra heating oc n° gives the best agree- 
ment for higher z values. The interpretation for 
the z heights < 800 pc have to be handled with 
care as this region is contaminated by radiation 
from the midplane and therefore the part respon- 
sible for the DIG emission is difficult to estimate. 
Magnetic reconnection (e.g. Birk et al. (1998) ) 
or heating by cosmic rays through linear Landau- 
damping (Lerche & Schlickeiser (2001) ) are possi- 
ble processes producing a heating independent of 
density. As photoelectric heating from dust grains 
is a n 1 , the data suggest that there is more dust 
present at z-heights below 800 pc than present in 
the models if this mechanism is responsible for the 
elevated temperature. 

In Fig. 7 the line ratios [Nil] /Ha and [SII] /Ha 
are plotted against each other. Values of [Nil] /Ha 
greater than 1, which cannot be explained by 
classical HH-region calculations, are reached with 
models using an extra heating source. Also in 
this case the extra heating independent of den- 
sity is able to produce higher [Nil] /Ha ratios than 
the other therefore matching better with the data. 
The limits of Haffner99 for two constant S + /S ra- 
tios (0.5 and 0.25) are also given. The ioniza- 
tion fractions of the models (see Fig. 5) are within 
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these two limits for z > 0.6 kpc. Together with the 
temperature plot (Fig. 4) the models match in all 
these direct and derived quantities with the data 
and agrees with the estimates of Haffncr99 and 
Reynolds et al. (1999). The application to other 
galaxies shows that this diagram is also a valuable 
diagnostic for the chemical evolution (Elwert et 
al. (2004)). 

4. Summary 

We have shown that the observed trend of the 
line ratios [Nil] /Ha and [SII]/Ha above the Galac- 
tic plane can successfully be explained by pho- 
toionization models including extra heating and 
considering the line of sight geometry. The ob- 
served values need an extra heating source which 
strength lies at the lower end of the predicted val- 
ues of Reynolds et al. (1999) due to our models. At 
high z heights (z > 800 pc) an extra heating inde- 
pendent of density gives the best agreement with 
the data, whereas for smaller z heights an extra 
heating term oc n 1 gives better results concerning 
the temperature. There is an intrinsically increase 
in the line rations and the gas temperature due to 
the progressive hardening of the radiation. The 
extra heating terms are enhancing this trend and 
elevate the line ratios to the observed values. It is 
important to incorporate the observing geometry 
to the Perseus Arm into the models when compar- 
ing with the data. A discussion concerning models 
of observed line ratios in edge-on galaxies is given 
in Elwert (2003) and Elwert & Dettmar (2005). 

This work was supported by DFG through SFB 
591 and through Deutsches Zentrum fur Luft- 
und Raumfahrt through grant 50 OR 9707. TE 
wants to thank Kenneth Wood and Ron Reynolds 
for helpful comments and enlightening discussions 
while writing the paper. We also want to thank 
the anonymous referee for making many very use- 
ful suggestions and comments which helped to im- 
prove the publication. 

REFERENCES 

Aller, L.H. 1984, Physics of Thermal Gaseous Neb- 
ulae (Dordrecht: Reidel) 

Birk, G.T., Lesch, H., Neukirch, T. 1998, MNRAS 
296, 165 



Bland-Hawthorn, J., Freeman, K.C., Quinn, P.J. 
1997, ApJ 490, 143 

Ciardi, B., Bianchi, S., Ferrara, A. 2002, MNRAS 
331, 463 

Collins, J.A., Rand, R.J. 2001, ApJ 551, 57 

Domgorgen, H., Mathis, J.S. 1994, ApJ 428, 647 

Dove, J.B., Shull, J.M. 1994, ApJ 423, 196 

Elwert, T. 2003, Ph.D. Thesis, Ruhr-Universitat 
Bochum 

Elwert, T., Dettmar, R.-J. 2005, ASPC 331, 203, 
editor: Braun, R. 

Elwert, T., Dettmar, R.-J., Tiillmann, R. 2003, 
A AS, 203, #111.05 

Ferland, G.J. 2002, Hazy, A Brief Introduction to 
CLOUDY, Dept. Phys. Astron. Internal Rep. 
96.00 (Lexington: Univ. Kentucky) 

Ferland, G.J. 2000, ASP Conference Proc. Vol. 216 

Ferland, G.J., Korista, K.T., Verner, D.A., Fergu- 
son, J.W., Kingdon, J.B., Verner, E.M. 1998, 
PASP 110, 761 

Ferland, G.J., Korista, K.T., Verner, D.A., Fer- 
guson, J.W., Kingdon J.B., Verner, E.M. 1998, 
PASP 110, 761 

Ferland, G.J., Mushotzky, R.F. 1984, ApJ 286, 42 

Haffncr, L.M., Reynolds, R.J., Madscn, G.J., 
Tufte, S.L., Jachning, K.P., Percival, J. P., 
Hausen, N.R. 2003, ApJS 149, 405 

Haffner, L.M., Reynolds, R.J., Tufte, S.L. 1999, 
ApJ 523, 223 

Hoopcs, C.G., Walterbos, R.A.M. 2003, ApJ 586, 
902 

Lanzafame, A.C, Tully, J. A., Bcrrington, K.A., 
Dufton, P.L., Byrne, P.B., Burgess, A. 1993, 
MNRAS 264, 402 

Lerche, L, Schlickeiser, R. 2001, A&A 366, 1008 

Miller, S.T., Vcillcux, S. 2003, ApJ 592, 79 

Miller, WW. Ill, Cox, D.P. 1993, ApJ 417, 579 

Mathis, J.S. 2000, ApJ 544, 347 



5 



Otte, B., Reynolds, R.J., Gallagher III J.S., Fer- 
guson, A.M.N. 2001, ApJ 560, 207 

Pauldrach, A.W.A., Hoffmann, T.L., Lennon, M. 
2001, A&A 375, 161 

Reynolds, R.J. 1993, in AIP Conf. Proc. 278, Back 
to the Galaxy, ed. S.S. Holt & F. Verter, 156 

Reynolds, R.J. 1990, ApJ 349, L17 

Reynolds, R.J., Haffner, L.M., Tufte S.L. 1999, 
ApJ 525, L21 

Sembach, K.R., Howk, J.C., Ryans, R.S.I. , and 
Kecnan, F.P. 2000, ApJ 528, 310 

Stafford, R.P., Bell, K.L., Hibbert, A., Wijcsun- 
dera, W.P. 1994, MNRAS, 268, 816 

Tullmann, R., Dcttmar, R.-J. 2000, A&A 362, 119 

Wood, K., Mathis, J.S. MNRAS 353, 1126 

Wood, K., Mathis, J.S., Ercolano, B. 2004, MN- 
RAS 348, 1337 

Wood, K., Loeb, A. 2000, ApJ 545, 86 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



6 





T 1 1 1 


t — ■ — i — ■ — i — ■ — i — ■ — i — < — 

fit of Haffner99 

standard model 

extra heating ~n° 1 

extra heating ~n 1 

WHAM data 


■ \ 






0^ 


K 






\ 
\ 

nig 














\ 






\ 



200 400 600 800 1000 1200 1400 1600 1800 2000 
vertical height z in pc 



Fig. 1. — Ha intensity of the Perseus Arm as 
observed by WHAM, overplotted is the fit of 
Haffner99 and the photoionization models with 
and without extra heating. 
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Fig. 2. — Sketch of the line of sight geometry and 
the beam smearing, which have to be considered 
when comparing the models with the data. De- 
riving temperatures from the data without taking 
this into account lead to temperatures that are on 
average 250 K higher. 
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Fig. 3. — The effect of the beam smearing and line 
of sight geometry on [Nil] /Ha specifically. The 
enhancement of this line ratio due to this effect 
is at most 15% and leads to higher temperatures 
derived from the observations if the effects is not 
taken into account. 
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Fig. 4. — Temperature structure of the DIG, es- 
timates from the observations with atomic data 
used by Stafford et al. (1994). The line of sight ge- 
ometry and beam smearing lead to temperatures 
on average 250 K lower, whereas the consideration 
of the ionization structure of iV + elevates the de- 
rived temperatures again. The figure shows that 
a combination of the two different heating terms 
gives the best agreement with the data. 
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Fig. 5. - Ionization fractions of N + and S + depen- 
dent on z, the assumption of Haffner99: N + /N = 
H + /H — ► 1 is approximately met for heights above 
1300 pc. 
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Fig. 6.— Observed [Nil] /Ha and [SII]/Ha in the 
Perseus Arm, showing the characteristic rise. The 
models with and without an extra heating source 
are able to explain this trend, the explicit values 
are better reproduced with an extra heating. 
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Fig. 7.— [Nil] /Ha - [SII]/Ha in the Perseus Arm 
as observed by WHAM, plotted are the models, 
the data, and the limits given by Haffhcr99 for 
two constant S + /S ratios. 
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